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ABSTRACT 


The generation of ambient noise by physical processes dependent on shearing of the 
boundary laver under sea ice is investigated. Special attention is paid to the identifica- 
tion of individual noise-generating mechanisms and the assessment of their relative in- 
portance. 

Recent studies of Arctic ambient noise are reviewed with specific reference to results 
showing particularly good or poor correlation between ambient noise levels and ice 
movement or relative current. Potential noise-generating mechanisms are described and 
categorized according to their small-scale driving forces and expected noise character- 
istics. More detailed quasi-objective investigations are then used to establish the relative 
importance of each mechanism as a contributor to the overall under-ice noise spectrum. 

Flow, Mechanical mechanisms, involving ice sheet fracture as a result of wind and 
current-induced bending moments, are found to be unlikely contributors. Conversely, 
processes in which ice fragments in current-driven motion under the ice interact to cause 
bumping and grinding noises. appear to be of probable importance. Turbulent pressure 
fluctuations in the boundarv laver under sea ice are shown to be of significance at low 
meauencies om a local scale. The role of resonant cavities im the under-surface of the ice 
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I. INTRODUCTION. 


Recent strategic thinking has placed great emphasis on anti-submarine warfare in 
the ice-covered waters of the Arctic Ocean (Young, 1986). This emphasis is the chief 
driving force behind ever increasing efforts to fully comprehend and model underwater 
acoustics in this region. A significant proportion of this effort has been directed towards 
the study of ambient noise under the ice and in the marginal ice zone. The importance 
of being able to predict ambient noise stems from its variability and its influence on 
passive sonar performance. 

The effectiveness of any form of passive sonar, be it shipborne, towed, air dropped 
or mounted on the sea bed, 1s described by the passive sonar equation. From Unick 
(1983) 


eee — NL — Di+ DI 


where 
SL is the source level: a measure of the sound emitted by the target 


TL is the transmission loss: the amount of attenuation experienced along the acoustic 
atin from the target to the receiver 


NL is the noise level: the level of background and self noise obscuring the signal at 
me TeECclVer 


DI is the directivity index: a measure of the directivity of the sonar equipment 

DT is the detection threshold: the minimum signal level that can be detected by an 

operator using the sonar equipment. 
It is the ability to accurately measure or predict the magnitude of the terms in this 
equation which permits the effective determination of passive sonar performance. A firm 
knowledge of the parameters which govern the passive sonar equation therefore allows 
an operator to plan both offensive and defensive operations with an invaluable degree 
of foresight. An example might be the ability to forecast a time period, or geographic 
location, offering poor passive sonar conditions in which a submarine might safely con- 
duct noisy evolutions while avoiding detection by hostile forces. Conversely, a unit en- 
gaged in offensive operations against submarine targets under the ice would benefit from 


an ability to predict a good acoustic environment. 


The magnitude of the source level term is dependent on the sound energy emitted 
by the object of the passive detection effort. It is therefore governed by engineering. 
naval architectural and tactical considerations not directly related to the environment. 
Transmission loss depends entirely on environmental factors, being a function largely 
of sea water and sea bed characteristics. Temperature and salinity gradients. both hori- 
zontal and vertical, are particularly important. Although this term is undoubtably im- 
portant in fixing effective sonar ranges, variability over moderate distances and time 
periods is relatively small in Arctic waters. This means that good climatological. or better 
still, recently measured, profiles of temperature and salinity are usually adequate for so- 
nar range prediction work over a broad area and for considerable time periods. 

The noise level term is generally of the most crucial importance in the overall sonar 
performance prediction effort in Arctic waters. It depends to a large degree on the en- 
vironment and is known to be highly variable both spatially and temporally. Noise level 
incorporates the effects of self noise, created by the hydrophone itself or by the platform 
on which it is mounted. and ambient noise from the surrounding water mass. Ambient 
noise from the environment is a variable factor depending on many environmental forc- 
ing functions. Its accurate prediction is crucial to the successful assessment of passive 
sonar performance. 

The last two terms of the equation, directivity index and detection threshold, are 
functions of equipment design. operating technique and in the case of the latter, opera- 
TOR Chiectiveness. 

Thus it can be seen that ambient noise prediction is an essential part of planning 
for Arctic operations. An additional utility can also be found mn weapons, sensor. and 
underwater communications svstems design. Where a grasp of typical ambient noise lev- 
els is important. 

Considerable effort has been directed towards the modelling and prediction of am- 
bient noise under ice. The problem has been tackled in a variety of Ways including em- 
pirical approaches, attempts to predict noise levels through the summation of emissions 
from sources and by consideration of energy budgets (Oard, 1987). Pritchard (1988) 
emphasised the need to investigate and understand the physics behind noise-generating 
processes on a local scale before modelling them and summing their effects at a partic- 
ular location. The overall predicted noise Would then be the sum of the noise intensities 
from many locations, generated by a variety of mechanisms and subject to transmission 
losses depending on receiver depth, source depth, range to the source and frequency. The 


resulting Arctic ambient noise model would be driven by outputs from existing 


tJ 


meteorological and ice dynamics models. A categorization scheme for nolse-generating 
mechanisms was proposed, by Pritchard, based on forcing factors. Mechanisms impor- 
tant under compact ice conditions are separated into three areas 

1. ridging effects which are driven by large scale ice stress and deformation 

2. micro-cracking effects which are driven by local ice stress 


3. effects dependent on the shearing of the boundary laver of the ocean under the ice. 


It is the purpose of this study to identifv and investigate any mechanisms which mav 
fall into the last category. Such potential noise-generating mechanisms are catalogued 
as exhaustively as possible and categorized according to their small scale driving forces 
and expected noise characteristics. Investigations are then conducted with a view to 
establishing their importance, if any, in the overall generation of ambient noise under 
ice. 

Manv of the physical processes involved are highly complex. and some are poorly 
understood or the subject of current research. For this reason the investigations are to 
some degree subjective. An attempt has been made to model the physics in as elementary 
a fashion as possible and generally to draw inferences from ‘order of magnitude’ results. 

Chapter II 1s a review of recent Arctic ambient noise research which considers the 
effects of boundary laver shear. The purpose of this brief review is twofold. Firstly to 
explore the fundamental concept that noise 1s caused by shearing of the boundarv laver 
under ice, and secondly, to identifv individual mechanisms which have been suggested 
as being responsible for noise. Evidence 1s sought to establish whether shearing in the 
oceanic boundary layer, in the form of ice movement through the water or water move- 
ment under stationary ice, 1s important as a forcing function for nolse generation. 

Chapter III outlines potential noise generating mechanisms and a proposed cat- 
egorization scheme. The individual mechanismis are then investigated in Chapter 1V with 


conclusions and a summary presented in Chapter V. 


Il. REVIEW OF CURRENT LITERATURE. 


Previously published studies of under-ice ambient noise are fairly numerous and in- 
vestigate a broad range of different forcing functions. Forcing bv thermal effects were 
investigated by Ganton and Milne (1965) and Milne (1972) while forcing by pressure 
ridging and ice mechanics was studied by Pritchard (1984) and by Buck and Wilson 
(1986). Various other atmospheric forcing functions have also been studied by Payne 
(1964). Ganton and Milne (1965) and Greene and Buck (1978). 

Three recent papers, Makris and Dver (1986), Lewis and Denner (1988) and Waddell 
and Farmer (1988). consider shearing of the oceanic boundary laver in the form of ice 
movement or flow of water under ice as a forcing function. These papers are therefore 
reviewed briefly with the objective of highlighting results which tend to confirm or deny 
the importance of boundary laver shearing under ice as a forcing function for noise 
generation. 

Makris and Dver (1986) analvsed low frequency ambient noise measurements and 
comprehensive environmental data from the 1982 FRAM IV expedition to the central 
Arctic (83°.\V 20°E£). The octave band 10 - 20 Fiz was taken as being representativemam 
noise in the low frequency range and compared with temperature and various applied 
stress components. The underlying hvpothesis was that:- 

Noise is created by ice fracture mechanisms proportional to the state of stress in the 
ice as induced bv environmental loading. 

Cross-correlation coefficients between the low frequency noise records and shear 
stresses due to the wind and current were found to be particularly high, 0.84. in both 
cases for 9.9-day records and 0.71 and 0.74, respectively, for 23./-dav records. {heveor 
responding coefficient for tensile stress due to cooling was 0.15. Composite measure- 
nents. namely ice stress on the ice sheet’s vertical section and the stress moment acting 
about the ice sheet’s central horizontal plane, were also highly correlated with ambient 
noise, correlation coefficients being 0.81 and 0.87, respectively, for 9.9 day-records. 

The general conclusion of the work was that low frequency pack ice noise correlates 
best with the moment due to opposing wind and current stresses acting on the ice and 


Worst with air temperature. 


Lewis and Denner (1988) studied the correlation between noise at 10, 32 and 1000 

Hz and a number of parameters describing temperature and ice kinematics. The data 
analysed were obtained during the AIDJEX project in the Beaufort Sea during 
1975-1976 and was exceptionally comprehensive. High qualitv positional data were 
available from some 40 satellite-tracked drifters which allowed the accurate determ- 
nation of ice translational speed and differential kinematic parameters. Correlation co- 
efficients were calculated for ambient noise levels at the above frequencies and the 
following parameters: 

eine translation rate of the ice 

= proportional to the kinetic energy of the ice 

eerie divergence Of the ice 

[D| 

fPecie vorticity (rotation rate) of thie ice 

IV 

Sethe normal deformation rate of the ice 

5. the shear deformation rate of the ice 

ieee temperature difference between the air and the ice 

and 


(GV? + S?)'?, the total deformation of the ice. 


The results of these calculations are reproduced as Figure 1. The correlation ob- 
tained throughout between ambient noise and LU, the simple translation rate of the ice 
measured by satellite tracking, 1s quite striking. In virtually all cases the correlation was 
the maximum observed for any of the parameters above. U correlated particularly well 
with 32 Fz noise in the summer (August 1975). This was attributed to the ice “rushing 
through the water”. Another observation was that. in general, LU correlated better than 
U2, except at 1000 Fz. This suggests that pure ice motion 1s more important in produc- 
ing noise than the kinetic energy of the floes, except at higher frequencies where collision 


mens are probably more significant. 
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Figure 1. Linear correlation coefficients between ambient noise and ice kinematic 


parameters (from Lewis and Denner, 1988). 


Overall this work seems to emphasize the importance of ice movement in the gen- 
eration of ambient noise. Although no particular mechanisms are mentioned, the role 


of shearing in the oceanic boundary layer appears to be of probable significance. 


Waddell and Farmer (1988) describe the noise resulting from the break-up of land 
fast ice in Amundsen Gulf in the Canadian Archipelago. Although the study was in fact 
motivated by the possibility that ice conditions could be monitored acoustically, the 
conclusions are of some interest here. 

Data collected over a 127-day period from April to August 1986 included ambient 
noise over the frequency range 300 Hz to 14.5 kHz, current measurements and 
meteorological observations. Currents in the area are strongly influenced by the tides. 
This 1s shown in Figure 2 which also depicts time series of longwave infra-red, ambient 


noise and air temperature for a 10-day period in June. 
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Figure 2. Time series of ambient noise and various environmental parameters 
(Waddell and Farmer, 1988). 
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Initially, when the ice was consolidated and land fast, there was no obvious corre- 
lation between the current meter observations and ambient noise. A more marked cor- 
respondence was noted between noise and atmospheric factors. A strong northerly 
current, developing in July, however, coincided with a sharp increase in noise levels as 


the ice near the recording instruments began to break up (see Figure 3). __ 
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Figure 3. Time series showing the increase in ambient noise levels with increased 


current and ice breakup (Waddell and Farmer, 1988). 


These results seem to indicate that when ice 1s landfast and static, a moderate cur- 
rent under it does not create discernable noise in the range 300Hz to 14.5 kIIz. When 
the ice cover is in motion, however, significant correlations can be found between cur- 
rent speeds and high frequency noise. This is probably be due to collisions and cracking 
during ice breakup. 

AS a general summary, the studies done by Makris and Dyer (1986) support the idea 
that low frequency (10-20 Hz) noise results from opposing wind and current stresses 
acting On the ice. The results of Lewis and Denner (1988) show strong correlation be- 
tween noise and ice translational speed which 1s proportional to the relative current un- 
der the ice. Waddell and Farmer showed that little noise results from relative current 


under stationary ice. 


WI. POTENTIAL NOISE GENERATING MECHANISNIS. 


Ambient noise due to shearing of the boundarv laver under sea ice could be caused 
by a number of different mechanisms. Little or nothing has been written about the de- 
tatled processes involved and in some cases the phvsical processes are sufficiently com- 
plex to prohibit meaningful analytical investigation. 

A number of potential noise-generating mechanisms are put forward here. The list 
is not intended to be exhaustive but most mechanisms which appear to have any possi- 
bility of being relevant to real world conditions are included. Some mechanisms appear 
intuitively more important than others while some may seem a little obscure. All are are 
included in this chapter and will be investigated to a greater or lesser degree in Chapter 
1M 


A. A CATEGORIZATION SCHEME. 

Two broad categories of potential noise generating mechanisms have been identi- 
fied. Firstly. Flow) Mechanical mechanisms which involve mechanical interaction between 
the sea and the ice cover and produce sound through cracking or impactive events. 
Secondly, Flow mechanisms which are purely a function of fluid flow in the presence of 
the solid ice boundary. 

The potential mechanisms are therefore divided and ordered in this study according 
to these criteria. Figure 4 summarizes the mechanisms that have been considered and 
their relative positions in this scheme of organisation. Figure 5 1s a schematic summiarv 


showing the partition of Flow’ Mechanical and Flow mechanisms. 
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Figure 4. An overview of potential noise-producing mechanisms. 
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Figure 5. A schematic summary showing the partition of Flow/Mechanical and Flow 


mechanisms. 


B. FLOW/MECHANICAL MECHANISMS. 
1. Cracking Noise Due to Wind/Current-Induced Moments. 

The concept of opposing wind and current stresses inducing moments in the 
Arctic ice cover was originated by Makris and Dver (1986). Although no detailed de- 
scription of the mechanism itself was presented, the high correlation obtained in the 
FRAM IV data, between noise and the stress moment acting about the ice sheet’s hori- 
zontal plane, indicates the potential importance of this mechanism. 

In an attempt to investigate the detailed mechanics of the processes involved, 
the wind and current stresses are taken to have most effect on features with some ver- 
tical extent. Pressure ridges or ice hummocks present obvious points where winds and 
currents could act to set up bending moments in the ice sheet which might possibly lead 
to ice fracture and therefore noise. The magnitudes of such bending moments are a 
function not only of the relative velocities of the Wind and under-ice current, but also 
of the size and shape of the ridge and of the thickness and general nature of the ice sheet 
itself. 

Maximum opposing wind and current velocities are most likely to occur when 
the ice 1s in a low-stress state with some degree of Wind driven motion. Ice in a static 
State is less hkely to be subjected to strong induced moments as the geostrophic and tidal 
currents in the open Arctic tend to be relatively weak. Wand and current stresses will 
also be greater in areas of highly deformed ice cover as drag there 1s more pronounced. 
Fracture as a result of induced bending moments may not neccessarily be more likely. 
however, as significant deformation is often coincident with greater ice thickness, lower 
salinity and therefore greater strength. 

2. Ice Fragment Bumping. 

The importance of ice collisions in the production of ambient noise has been 
recognized in several works (Diachok and Winokur .1974; Shepard, 1979) and is an es- 
tablished mechanism when considered on a large scale, 1.e., floe-floe interactions. The 
presence and movement of significant amounts of ice debris under sheet ice 1s also doc- 
umented (Zubov, 1943, Buck and Wilson, 1986). A potential noise generation process 
therefore exists in any pack ice covered region where deformation events create quanti- 
ties of loose fragmented material which can be moved around by a relative current. 
Noise from this source is likely to be aurallv detectable and characterized by bumping 


Or grinding sounds. 


Ihe) 


C. FLOW MECHANISMS. 
1. Noise From Turbulence. 

Turbulence is a recognized source of low frequency noise in the open ocean 
(Wenz, 1962; 1972). Although the processes by which it is radiated are exremely ineffi- 
cient and its far field effects are minimal, the direct effect on hydrophones within tur- 
bulent flows can be significant. This is because the pressure fluctuations impinge directly 
on the active surfaces of hydrophones which react by generating an electrical signal. 
Thus turbulent pressure fluctuations can cause a type of self-noise whose effects mav 
be significant at very low frequencies. 

In the Arctic the situation 1s complicated by the presence of ice cover. The ice 
acts as a boundary of variable roughness past which water flows. In addition. the ice 
sheet itself mav play a role in the radiation of turbulence noise to ranges bevond the 
actual source region. Turbulent self-noise has been recognised as a possible contaminant 
of Arctic ambient noise data by Makris and Dyer (1986) and is included as a noise- 
generating mechanism in its own right. 

Flow-induced noise and flow-excited noise are described briefly in Chapter 1V 
as potential mechanisms by which the effects of turbulent pressure fluctuations are ra- 
diated to the far field. 

2. Resonant Cavities. 

The under-surface of floating ice in Arctic regions is, by its Verv nature, rough, 
deformed and pitted with cavities caused by melting, fracture or brine rejection. Such 
holes or cavities present a potential noise source when water is in motion across their 
openings. Flow mav be the result of ice motion or currents under static ice and could 
excite resonant responses in these cavities. 

No known reference has been made to the possibility that resonant cavities 
might be responsible for under-ice noise but similar above-water features are certainlv 
capable of producing noise in strong winds. Noise from hull cavity resonance 1s also well 
known in naval architecture (Blake, 1986). 

The parameters affecting the frequency and intensity of sound from resonant 


cavities Would be cavity size, shape and water flow speed. 
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IV. AN INVESTIGATION OF POTENTIAL NOISE MECHANISMNIS. 


A. FLOW/MECHANICAL MECHANISMS. 
1. Wuind/Current-Induced Moments at Pressure Ridges. 

The action of a current relative to an ice field on submerged pressure ridge keels 
and the action of the wind on their above water sails gives mse to a potential 
flow. mechanical noise-generation mechanism. The moment set up by these forces may 
be enough to cause cracking and cracking noise in the ice sheet close to the ridge fea- 
UTES. 

The moment j/, due to water pressure on the ridge Keels is likely to be rein- 
forced by wind action on corresponding sail formations causing a moment JJ, (see 
Figure 6). This is particularly so when internal stresses within the local ice field are low 
and motion is largely due to wind-induced free drift. The applied moment J/, is therefore 
a function of the wind velocity J), . the relative current velocity 1”, and the pressure ridge 
Peemetis. if, is the sum of 3/, the moment caused by the relative current on the keel 
Aeerty, the moment on the sail due to the wind. 

This functional relationship is investigated below by assuming an idealized 
Meesstire ridge model after Wright et al. (1978) and Paquette and Bourke (1988). The 
geometrv of this model is illustrated in Figure 7. 

It is assumed that 


pete SvStem is in localisostatic equilibrium when V, = = 0 


2. the ridge is svmmetrical about point O as shown in Figure 7, with dimensions 
governed by local isostatic balance 
3. the density of the ice is homogeneous throughout with an equal percentage of void 


in the sail and keel structures 
4. the surface of the ice sheet. sail and keel are perfectly smooth 


>. the wind velocity v,(z) and relative current velocity v,(z) varv linearly with height 
and depth as shown in Figure §&. 





Figure 6. The action of the wind and relative current on a pressure ridge feature. 





Figure 7. The idealized pressure ridge model. 





Figure 8. The assumed linear wind/current profile. 


a. Ridge Model Geometry. 


From the assumption of isostatic equilibrium 
Vals W.p;= (tg, + Yahg(I, — WM Ow— pd (1) 


where 
p, \S the water density 


p, is the ice density 
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from Figure / 


ie ? 
— tance, (- 
and 
—_ 3) 
we tan 2. ( 
SO 
(/ aa 4 
k tan Cp cs 2) ( ) 
Substituting (2) and (3) into (1) and rearranging 
he ~; oy Ps tana, 
hy ee 5 7 hy, tan je 
| a ipo 
Solving this quadratic in x, ) 
ly 
: 
~; tan o. : 
I+; 0+ —-—— a Ue 
We Py pi tan oj 
oe =e ee ee ae (S) 


Ny 7 p; 
Pw — Pi 


Now taking tvpical values for p, and p, as 920 kg m3 and 1026 kg m-, respectively, and 
feo) oc. = 20° (from Wright et al., 97S) 


0.4) 6 
P=0 (6) 
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also from Figure 7 


‘g 


Ay 
hy, = lr + hy + a tan Ay 


Or 


tan Op 


eet 
k ag me baler: 


+ ty (7) 


Since the ice sheet itself can also be assumed isostatically balanced 


50, = (Py — Pi) (8) 
and since 

f=1,+ t, 

n= 1th = 0.91 (9) 
and 

pm OE sim (10) 


Substituting (9) into (7) 


tan % pp Pp; 


ly, = hh, +A, —— 
Ke K zi 5 tan oie Pw 


Further substituting numerical values for a, «, and h, from (6) 


= ayiae th. (11) 
Or 
fh ae eas + fp, (12) 


From the geometry of the sail 
Pm fet, (13) 


and from (10) 


; (Qa pi) 
in. = ia ae ie 
SO 
h’~h,+O0.12 (14) 


Now that functional relationships have been established between the vari- 
ous ridge dimensions, an equation for \f,, the moment caused by wind action on the sail 


about point O can be found (see Figure 9). 
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Figure 9. Wind action on the ridge sail. 


a) 


b. The Moment Caused by Wind Action on the Sail. 
Using concepts from Zubov (1943) and the theorv for a fluid jet acting on 
a fiat plate (Hannah & Hillier. 1970) the wind velocity v, at height 2 above sea levemeam 
be resolved into two components: v, the component parallel to the surface of the ridge 


sail and r, normal to this surface. 


aco aee (15) 
b= hy, Sin o, (16) 


At the windward surface of the sail », causes a dynamic force. v, has no effect if the sur- 
face of the sail is perfectly smooth. The resulting force P is equal to the loss of momen- 


tum to the sail per second. So 


=) 


~ 


P=nv, sing, kgms 
Where 

21is the mass of fluid (air) incident on the sail at height z. per second. 
iihen 

7 ae -2 ‘ 

P =p, Sit a (gees (1% 
where 

pais tne air density 


The moment 7,(z) about point O is given by the product of the length of the 


lever arm @ and the component P,. of P. normal to @ (see Figure 10). 
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Figure 10. Moment about O, the point of symmetry, due to the wind. 


From Figure 10 


P,=P, cos B 
and 
B=n/2-—O0+0, 
SO 


P,= P sin(@ =— a.) 
and from (17) 


D. . : 
P, = pqv, sin a, sin(@ — a,) 


23 


(18) 


Now sin(@ — «,) can be found in terms of f,’, ?, a, and z by considering triangle OAB, (see 


Figure 11). 





Figure 11. Triangle OAB. 


From Figure 11 


sin(@ — o,) = sin @ cos a, ~ cos @ sina, 


or 
X COS a, —Z Silla, 
sin(@ — o,) = aaa eam 
where 
, 
a ee 
tan o, 
SO 
nee 
COS G2 
Ste 6) =e 


¢ sin &, 
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(19) 


Substituting (19) into (18) 


2 2 
_ Patifs’ COS"a, ~ =) 


pet (20) 
and 
m(z)=¢P,= pra(h,! COs*a, = 7) (21) 


If the wind shear ts linear, as outlined in the assumptions above 


Vz — 6) 
= (22) 
Sy 
SO substituting (22) into (21) 


<3 
Onin 





mz) = 2 = ete COs’ ¢, =) (23) 


A. 


The total moment ./, on the sail about point can now be found bv inte- 


grating with respect to z between ¢, and A,’ 


= ee 
os ae hee | 
Mo=— | (2 -1,)'(hy cos*e, — 2d: (2-4) 
Ay 
ts 
OF. 
1’ 
(Ue 2 5 2 : Ps, 5 
Me i : 2 (h, “(4 COS‘ or, ~~ 3) =e h, 4G | cosa) oF h, *12(6 = 12 cosa.) — dh, ic cosa, + 1") 
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From (13) 


h’ =h, +t, 


SO 
Pal 2 2 
M,= 5 2 (he (4 cos"a, —3) + hy,(4 cos’a, — 4)) 
or 
One ye: 25 
We = pV AteG. (25) 


where C, is a form coefficient dependent on the geometry of the ridge. 


/ A I, \ 
C,= as (2 COC 1 (4 cos’, — 2-9) (26) 
Lio 20 cecig =—Ome 
Gee (4-456 x 107° (27) 
He 


Figure 12 shows C, plotted against h, for various values of t. Note that C, ~ 0.04 for sail 


heights greater than one metre and 1s virtually independent of 7. If o,, the slope of the 


Ne 


sail. 1s allowed to vary, an interesting effect 1s noticeable (see Figure 13), 9 item 


+ 


o,2 28° , C, reverses sign, 1.c., the result is a moment in the opposite directions 
could be relevant when considering young, relatively unweathered ridge features where 


slopes mav be greater than the 20° assumed above 
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Figure 12. The form coefficient for a pressure ridge sail as a function of sail height 


and ice sheet thickness. 





Figure 13. The form coefficient for a pressure ridge sail as a function of sail slope. 


2/ 


c. The Moment Caused by Current Action on the Keel. 
The moment 4/, caused by the relative current on the keel can be found in 


a similar way to that found for the sail. The relevant geometry is outlined in 


Figure 14. 





Figure 14. Geometry for calculating the moment caused by relative current on a 


pressure ridge keel. 
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By comparison with (21) 
ee? 
i oe a ae = Pw (/, COs Op — Zz) 


Linear current shear 1s assumed, so 
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ee 
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p @ 2 
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af 
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ns 


Meo = 30° and ¢, ~ 0.97 
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My = = (2h; — 0.92,0) 





or 


yo7.2 
My = om oly Cr 


were ©, 1s a form coefficient for the keel given by 


C= (163-75 ) x 107 


Ue 


Figure 15 shows C, plotted against h, for various values of t. There 1s more dependence 
on ¢ than in the case of C, and it should be noted that C,=O°13 for 7) — iit 
comparing Figures 12 and 15. C, is generally about four times bigger than C,. When 
o, 18 increased, there is a point at about 44° where VJ, reverses sign as in the case of WA 
at greater slopes (see Figure 16). Again this may be relevant when considering relatively 


new, unweathered, ridges. 
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Figure 15. The form coefficient for a pressure ridge keel as a function of sail height 


and ice sheet thickness. 
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Figure 16. The form coefficient for a pressure ridge keel as a function of keel slope. 
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The total moment J/, on the ridge structure about point O is given by the 


sum of the moments caused by wind on the sail and relative current on the keel. 
M,= M,+ My, (37) 


Figure 17 shows J, with its constituent components Jf, and 37,, as given by (25) and 
(35), plotted against , with ¢ = 1.5 m. Typical values for p,, V,. p, and V, have been used, 
they are summarized below:- 

p,=1.3kgm°> 

p, = 1026 kg m3 

= 10m ss 


(0 is 


These wind and current velocity values were chosen to represent the upper limits of ob- 
served measurements from the hterature , i.e., a worst case has been chosen in order that 
the resulting moments will be the largest that could normally be expected to occur in 


hature. 


gth) 


& 
8 
re 
i= 
= 
3 
5 
ec 
eg 
2 
o 
3 


h (metres) 





Figure 17. Components of the applied moment at a pressure ridge on 1.5m thick ice 


(kKNin per m ridge length). 
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fyeure 18 is a plot of M, against ¢ and f. with V,=l10ms7? and V, =0.2ms"'. The 


weak dependence on ice sheet thickness is apparent as is the order of magnitude of MM, 


to be expected for these conditions , 1.e., Af, ~~ 1 KN m per m ridge length. 


t (metres) 


=a 
So 
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Vw = 10m/s, Ve = 0.2m/s 





The total moment applied by wind and current at a pressure ridge as a 


Figure 18. 
function of sail height and ice sheet thickness (KNim per m ridge length). 
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The following points arise from the preceding calculation of the moment 
due to wind and relative current action on a pressure ridge. 
1. Wf, is the dominant component of .5/,. 
2. ¢, the thickness of the ice shect has little effect on the absolute value of 47. 


3. For angles of slope o, and o, taken from Wright et al. (1978) 37, is anticlockwise 
when the wind blows from the left and the current acts from the right (as shown in 
Figure 6). If steeper stopes are assumed, .\f, could be reversed. Regardless of the 
sense of the resulting moment, however, the absolute value of \/, at the ridge will 
Bevo much the same orden, 


There are two considerations which have not been taken into account when 


calculating /, and \/, which might lead to increased moments at the ridge. 


1. The surfaces of both the sail and keel have been assumed smooth. If a more real- 
istic approach were to be taken and some roughness allowed, drag could be ex- 
pected parallel to the sail and keel surfaces. 


ty 


Reduced pressure in the flow downstream of both the sail and keel would also have 
some effect on Jf, and J,. This would contribute in a clockwise sense for ama 
blowing from the left and current from the right. 

Despite these possible inaccuracies and any others which might arise from 
the assumptions made, equations (25) to (28) and (35) to (37) allow a reasonable esti- 
mate of .\/, to be made for sensible values of ridge height and ice sheet thickness. 

2. The Effect Of An Applied Bending Moment On The Ice Sheet. 

In order for the moment J/,. described above, to cause ambient noise undemume 
ice, cracking must occur in the ice sheet. The horizontal distribution of bending moment 
Af{(x) on the ice sheet must therefore be investigated in order to estimate the maximum 
moment that could be expected under normal circumstances. Two extreme cases are 
exanuned to determine whether the maximum expected bending moment on the ice sheet 
is likely to be sufficiently large to cause cracking. The first case examines the ridge as 
a relatively new feature when it 1s in an unconsolidated state. That 1s, it is assumed that 
the ridge structure consists of loose rubble and ice blocks and therefore contributes no 
strength to the system. The second case, representing an older pressure ridge, occurs 
when the ice debris comprising the sail and keel has frozen into a consolidated mass 
which together with the ice sheet near the ridge can be considered rigid. Figureme 


compares, schematically, the two conditions outlined above. 
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Figure 19. Consolidated and unconsolidated pressure ridges. 
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a. The Unconsolidated Pressure Ridge. 

In this case the ridge ice sheet svstem is assumed to behave in the same 
Way as two adjoining. loaded, semi-infinite beams on an elastic foundation. The uncon- 
solidated mass of rubble and ice blocks forming the sail and keel interact with the wind 
and current as described previously to cause an applied moment A/, and also a distrib- 
uted loading on the ice sheet. Thev do not however contribute in any way to the strength 
of the sheet. In order to calculate the expected bending moments on the 1ce sheetsmae 
section to the right of the ridge axis is treated as a semi-infinite beam with static loadings 
contributed by the sail and keel as shown in Figure 20. The total loading can be thought 
of as being due to three triangularly distributed loadings. That is, a downward load due 
to the weight of the ice in the sail (triangle ABC in Figure 20) and an upward loading 
due to the buovancy of the keel (triangle DEH minus triangle FGH). So in the isostatic 


condition 


YAN) = Vase —¥pen t+ ¥FGH (38) 


Where vis the downward displacement of the ice sheet. or 


silx) = J (downy ~ Sup) (39) 
Similarly 

Mix) = Masco — Moen + Mren (40) 
or 

AV Conic eee or (41) 


M is the bending moment on the ice sheet (positive in the clockwise direction). 
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Figure 20. 
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Static loading on the ice sheet by an unconsolidated pressure ridge. 
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Equations for the displacement and bending moment on a semi-infinite 
beam subject to triangularly distributed loading have been developed from Hetenvi 


(1947) in the appendix. They are:- 

MX) = QoF lA, x) 0<x<7 Gm 
where F(/, ?, x) = (Cy,_,, + 42(2 — x) + BA, — D,,( Aye +2B,2)) 

H(X) = GG(4, €,x) x2? (43) 
whére G(/, 7, x) =(C,,..) + 2.47 — Oe ee 

M(x) = — qo H(A, @.x) O<x (44) 


where H(/,?.x)= (A ites = DiAy— Bb Ap ee) 
¢ 1s the base length of the triangular loading 
and g, 18 the loading at the free end of the beamiere aig —0 


see the appendix for definitions of 1/4 the characteristic length, k the foundation 
modulus and the functions 4,. B,. C, and D,. 


The net downward loading gyq,,,) 1s found by considering triangles ABC and 
PGi (See Picure 20): 


Focdowny = MsPiS + (My! — hy — HM Pw — PIS (45) 
and using (11) and (6) 

Jo(downy = Msg; + 1.6(0, — 3) 

Taking 9-= 920 Ke nic and 9p. = 1020 ken 

Geen = 0G 0 acon (46) 
From Figure 20 


Z Ses 
(down) " 2 ~ tan a, 





~ 2.8h, (47) 


38 


For the upward loading ¢,,,) (due to triangle DEH) 
Youp) = hry — pie 
Fo(up) = Soy os 


or, using the constant values given above 


- (8 
fey 4155h, kgm ss ~ 


fmeomrriure 20 





. 
a 2 tana, * ? = 1h, 


(48) 


(49) 


(50) 


Substituting for gq and ¢# in (42). (43) and (44), also taking A =(p, — p,)g = 1039 


mem -S~? 


ee) = 0.942 FA. Lpaowny X) 
roan) = 0.942 GUA, Caowny ©) 
en) = —459 HG, fgowny- ©) 
Vay) = 0.144 FA. Leap. ©) 
Yup) = 0.144 Gz, Zan. *) 


Mp) ) = —75 (A, ne x) 
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Oa yeaa) 


X= t(down) (2) 


Both @,, and &4...) are related to h, by (47) and (50) and / is a function of 7, so from (39) 


VAX) = 01942 Fedownyt le yi) — Volt hae ore) OSX < CO gowny ee 
r(x) = 0.942 Crees h.. x) ~ 0.144 oa ha ¢ (down) = ne a (58) 
De) = OFA NG oui ge) Ge eee xX > Cup heed 


and from (41) 


Mx) = —485 Hegowny(t. Agex) + 75 H 


tupy\ts he) O<x (60) 


down) 


Figures 21 and 22 show the displacements of ice sheets of various thick- 
nesses When subjected to the loading of unconsolidated ridges with sail heights of 1.5m 
and 3m. It can be seen that the concentrated mass of the sail causes downward dis- 
placement at the ridge axis. The nore distributed buoyancy force due to the keel camaes 
upward displacement awav from the ridge. This results in bending moment curves as 
shown by Figures 23 and 24. The value of the constant 4 was calculated using a tvpical 
value for E, the Young's Modulus for ice. E = 3 x 110° Pa, irommewellor an 

The absolute value of the bending moment due to the isostatic loading on 
the ice sheet reaches a maximum value JV/,,,, some 10 to 20 m from the ridge axisauiie 
at this pomt, with or without the additional apphed moment Jf, due to the wind and 
relative current, that cracking 1s most Itkely to occur. Figure 25 shows 34%). plgiias 
against A, and ¢. 

Having established estimates for Jf, the apphed moment and \,,, the 
maximum bending moment on the ice sheet due to isostatic loading, the two can be 
compared over a range of ridge heights and ice sheet thicknesses. Fron: Figure 1$ and 
Figure 25 it is apparent that ./, is of the order 0 to 1 KN m per m ridge length fommeae 
istic Maximum wind and current speeds. .M,,,,, however, 1s of the order 20 to) 2000RRas 
m per m ridge length, i.e., \/,,, 1s roughly two orders of magnitude greater than JV, for 
a given ridge size. 

Thus, it is concluded that if cracking occurs in the unconsolidated case, it 
will be due almost entirely to the isostatic loading on the tce sheet. The additional ap- 
plied moment JZ, while not being sufficient to cause cracking by itself, could possibly 
act as a trigger mechanism for cracking where a cracking event due to 1sostatic loading 


was likely to occur anyway. 


40 


40 50 60 
x (metres) 


Figure 21. Ice sheet displacement due to isostatic loading of an unconsolidated 





pressure ridge, height 1.5m. 
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Figure 22. Ice sheet displacement due to isostatic loading of an unconsolidated 





pressure ridge, height 3.0m. 
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Figure 23. Absolute bending moment on the ice sheet due to isostatic loading of an 


unconsolidated pressure ridge, height 1.5m. 
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Figure 24. Absolute bending moment on thie ice sheet due to isostatic loading of an 


unconsolidated pressure ridge, height 3.0m. 
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In order to establish the likelihood of cracking due to isostatic loading 
alone, Af, 18 compared to the critical bending moment AZ,,, at which cracking will oc- 
cur. A/,,,, 1s dependent on a, the tensile strength of the ice sheet, and its thickness t. A 
nominal value for the tensile strength of the ice sheet was taken, o, = 0.4 x 10° Pa from 
Mellor (1986). 


G,t 


Merit a 6 (6 1) 





Figure 26 shows the ratio A,,,,/M.,, plotted against A, and t. From this graph it can be 
seen that Af,,,/ALl 


crit 


~l where t~h,. When S7@,,,/1,,,21 cracking due to isostatic load- 
ing is likely to occur. It should be noted that as A, becomes much greater than ¢ the result 


is a greater likelihood of cracking. 





Figure 25. Maximum bending moment on the ice sleet as a function of sail height 


and ice sheet thickness (KN m per m ridge length). 
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Figure 26. The ratio between maximum bending moment and critical bending mo- 


meut for an unconsolidated ridge under isostatic loading as a function of 


sail height and ice sheet thickness. 
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b. The Consolidated Pressure Ridge. 

In this case a completely rigid svstem 1s envisaged. It is assumed that the 
ridge structure has undergone repeated freezing events and 1s no longer active in that 
ridge building processes have ceased. The material comprising the sail and keel has fro- 
Zen into a solid mass and contributes fully to the strength of the ice sheet. The solid ridge 
section 1s therefore treated as a rigid beam pivoting about point O (see Figure 27) while 
the adjoining ice sheets either side of the ridge are assumed to act as semi-infinite beams 
on elastic foundations. The boundary at section A-A, between the rigid portion (the sail, 
keel and an unknown length of ice sheet) and the rest of the ice sheet (which 1s consid- 
ered flexible), es a distance L from the ridge axis. Equations for the displacement and 
slope of a semi-infinite beam on an elastic foundation are taken from Hetenvi (1947) and 
can be found in the appendix. 

For simplicity, symmetry about the ridge axis is assumed such that 


iene distances 1, and 1, equal. 


tJ 


my eandey., the displacements of the ice sheet at distances L either side of the mdge 
axis, are equal and opposite. 


acne, tie Dendine moments on the ice sheet at distances L either side of the 
ridge axis, are equal and opposite. 


re 


4. QO, and Qa, the shear forces at distances L either side of the ridge axis, are equal and 
opposite. 


Thus 


Dy — —Lp = —[ 
Avi = —\'p =a 


M,=—Mp=M 


and 
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Figure 27. Geometry and general structure of the consolidated pressure ridge model. 
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Now summing the moments about point O 
LANo = M; ar My, ae Mp ome Oy aa LOp ae AY = 0 
or given the assumption of svmmetrv above 


M,+ 23M +2LQ =0 (62) 


AM, 1s the moment caused by the immersion of ice to the left of O and emmersion of ice 


foethe right. 


iL 
M,= | x*ktan@dx =2/3kL’ tané 
—L, 
Ol 
M, = 2/3KL°@ for small (63) 


From the appendix, the displacement of a semi-infinite beam subjected to a shear force 
wend moment ./ at its free end rs given by 
207 2g" 


os — 0: : 
‘ k ae k 





Cr 
See the appendix for definitions of 1// the characteristic length, k the foundation 
modulus and the functions 4, B, C,and D, . Since the section A-A, distance L from the 
meee axis, comcides with the begining of the flexible tce sheet and therefore the [ree end 


m@imune semi-infinite beam model, x =O and D, = C,.=1. Therefore 





Nie DI ae, | 
y= : Snaey (O — MA) (64) 
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The slope @ of the semi-infinite beam section is also given in terms of Q and .V/ in the 


appendix. 





20+ Vis 
k 


C= — Avert I: oe 


AGain at the frce ena <— 0 sowie —— aa 


20) aypi° 55? 
g-- 284 AME | BD ony oy (65) 
k k k 
Also, from Figure 27 
tan @ = ae or t=O, for small? (66) 


1S 


Now subsituting (63) and (65) into (62) 


= ECON eNO LLCO = ip 0 (67) 


and from (64). (65) and (66) 


eee) a 
~  (1+22) fo 
S067) becomes 


MCS GLY S4GLY 2402) Se 


uo { 4 


The effect of the applied moment J/, on \/, the bending moment at the interface region 


between the rigid pressure ridge section and the flexible ice sheet, is therefore given by 


(1+ (4L)) 
oe ae 7 
Ma Cao) + 4(AL)" + 4{AL) +2) 


Figures 28 and 29 show the ratio 1//\/, plotted for sail heights of 2 m and 4 m with 
various ice sheet thicknesses. 4 has been calculated using a Young’s Modulus of 
E=3x 10? Pa and wind and current speeds of 10 ms and 0.2ms~ have againsbeem 


assumed to represent a worst case. 
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Figure 28. The ratio of applied moment to the maximum bending moment on the ice 





Sheet for a consolidated pressure ridge, height 2m. 
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Figure 29, The ratio of applied moment to the maximum bending moment on the ice 


Sheet for a consolidated pressure ridge, height 4m. 
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i : 
Note that (69) is not valid for L< ar since the keel and/or the sail struc- 
tures significantly increase the rigidity there. It is assumed that cracking is unlikelv 
Within the ridge structure itself due to its increased thickness and strength. From 


SS 29 the maximum bending moment on the ice sheet itself appears to occur where 
k 
= 

The maximum bending moment ,,, 1s plotted for various values of sail 
height and ice sheet thickness in Figure 30. Note that A7Z,,, 1s of the order 0 to 0.5 KN 
m per metre of ridge length. Comparing this with A/,,,,, the bending moment required for 
cracking to occur as illustrated by Figure 31, the bending moment available is about 
three orders of magnitude smaller than that required for cracking to occur in the ice 
slieet. 


Thus, it 1s concluded that cracking noise is unlikely to occur due to 


wind/current induced moments applied to a consolidated pressure ridge. 
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Figure 30. The maximum bending moment on the ice sheet due to wind/current- 
induced moments at a consolidated pressure ridge (KN m per im ridge 


length). 
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Figure 31. The critical bending moment required for cracking to occur In ice as a 





function of thickness (KN m per m ridge length). 


3. Conclusions : Cracking Noise Due to Wind/Current-Induced Moments. 

From the results of a) and b) above it is apparent that the generation of noise 
due to these processes 1s unlikely to be of any major significance, if it occurs at all in 
nature. In the case of the young, unconsolidated, pressure ridge it appears that the 
bending moments on the ice sheet due to the isostatic loading of the sail and keel 
structures far outweigh any additional bending moments caused by the effects of wind 
and current on the ridge. Any cracking which occurs 1s likely to be due to these loadings. 
An additional contributor to local ice sheet bending might also be large scale ice stress. 
This has not been considered here as it would probably be small if the ice were in wind 
driven motion. There remains the remote possibility that the cracking of heavily loaded 
ice sheets could be triggered by the onset of wind/current-induced effects. It seems un- 
likely however that this could account for any significant correlation between wind, 
current and ambient noise in field measurements. The results of the investigation into 
the more mature consolidated pressure ridge reveals even less probability of cracking due 
to wind/current induced bending moments. The strength of the ice is such that it is very 
unlikely to yield to the relatively Weak bending moments so induced. This conclusion is 
Supported by the fact that significant noise levels have never been observed emanating 


from old inactive pressure ridge features. 
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4. Ice Fragment Bumping. 

This mechanism embraces all the bumping and grinding noises which would be 
expected to occur when ice fragments move about under an tce sheet. The probability 
that this mechanism 1s responsible for the creation of significant levels of ambient noise 
1s reinforced by the recorded observation of audible sounds from the ice surface. Zubov 
(1943) quotes Weyprecht as stating that 

Sometimes shift at the bottom [ of the ridge keel ] 1s audible, with the complete re- 
pose of the ice at the top. This occurs probably as a result of the movement of water 
under the ice field. The difference of the movements of the ice field and of the waten 
on which it lies, that is, the current of the water, is that power which levels the longen 
depth of the ice. 
Zubov (1945) also quotes Makarov on the existence and movement of ice blocks under 
te 1ce 
During the third Winter, a fissure was formed under the strip under the Fram. The 
fissure opened. lumps of significant size began to float out from below. “This 
shows”, savs Makarov, “that many lower lumps constantly travel. The current of the 


water and the movement of the ice change their direction so that the migrating ice 
block stops under certain conditions, under others it can move from the spot.” 


It is neccessary to add that the lumps which compose the under-ice part of the 
hummocks [ ridges ] are not onlv transformed but are also gradually destroved. 

Ice blocks or fragments are most likely to be the products of active ridging 
events. In the previous section the behaviour of young unconsohdated ridges was ex- 
amined quantitatively. The downward depression of matertal at the ridge axis was shown 
in Figures 21 and 22. This fractured ice appears to be thrust downward and trapped in 
buovant equilibrium under the ice sheet by the ridging processes. Buck and Wilson 
(1986) describe the forcing of ice blocks as large as four metres thick into the water 
column in their study of an active pressure ridge. 

It is reasonable then that given ice sheet motion relative to the water body be- 
low, the loose material will move horizontallv. Indeed, further observations bv Makarov, 
(Zubov, 1943) suggest that this is the process by which pressure ridge keels become 
vertically truncated and more widely spread in the horizontal than the corresponding sai 


formations. 


a 
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Given then that there are large ice fragments in motion along the rough under 
surface of an ice sheet. collisions, grinding and bumping are certain to occur. Noise from 
this source would probably be similar in nature to that observed as a result of floe-floe 
interactions in the marginal ice zone. Predominant frequencies are in the range 10 [1z to 
1 KHz according to Diachok and Winokur (1974). Certainly the impactive nature of the 
processes and the variety of fragment sizes is likely to yield noise spread over a wide 
frequency range. 

The intensity and perhaps the frequency range of ambient noise caused bv this 
mechanism would be dependent on 

mee ihe relative current Under the ice. 
Peele proferation of ice debris under thie ice sheet. 
3. The roughness of the under-ice surface. 


4. The thickness and density of the ice sheet. 


The current strength 1s of prime importance and 1s directly related to the translational 
eects Of the ice and or oceanic or tidal currents. The presence of debris and ice 
roughness are both a function of the deformation history of the ice field and are hkelvy 
to vary by region. Item 4 has been included as these factors are almost certain to have 
a bearing on the nature of the sound resulting from impactive events and its propagation 
into the sea. 

Modelling this mechanism 1s likelv to be a complex task if approached from a 
purely physical standpoint. A more profitable approach might be to seek empirical re- 
lationships between the factors outlined above and measured ambient noise levels. This 
could be backed up by some simple field experiments involving perhaps acoustic moni- 


toring of ice blocks pushed through bore holes. 


B. FLOW MECHANISMS. 

The potential noise creation mechanisms falling into this category are all driven ex- 
clusively by the interaction of the sea with the underside of the ice sheet floating at the 
surface. The factors important in their operation are therefore the flow velocity of the 
water relative to the ice, the characteristics of the under-ice surface including weaves 
shapes and sizes, the roughness, and to some extent the thickness, of the ice cover. 

It seems intuiuively obvious that some noise. no matter how small. must be 
produced when water flows under rough ice. The object here is to explore the mech- 
anisms Which may be responsible for this flow noise and assess their importance in the 
overall generation of ambient noise in an ice covered environment. 

Because of the complexity of the mechanisms outlined in this section, involving as 
they do, turbulent flow over largely random surfaces, the prediction of the noise fields 
created is somewhat crude. Some of the physical processes are currently being researched 
and many of the concepts are stil: only poorly understood (Ross, 1987). Where possible. 
estimates are made of likely sound characteristics in terms of frequency and intensity. 
These estimates are based on objective considerations but must be assumed accurate 
only to an order of magnitude. 

A typical composite Arctic ambient noise spectrum is presented in Figure 32. This 
serves to indicate the intensities required in order for a particular mechanism to be con- 
sidered important as a contributor. It 1s probable that flow noise, whether or not it is 
found to be a significant component of the total noise under pack-ice, would be entirely 
masked bv noise from other mechanisms at locations near the ice edge. For this reason. 
the measured spectrum chosen to represent the norm is from a central arctic pack-ice 
location. It 1s reproduced from Makris and Dyer (1986) who made measurements at 
So 20 Ba Apre ls 
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Figure 32. A typical under-ice noise spectrum (from Makris & Dyer, 1986). 


1. Noise From the Turbulent Boundary Layer. 
a. Boundary Layer Characteristics. 
The ratio of kinetic to viscous forces in a flow regime is described by the 
Reynolds number for that flow. 
ee L, | 

Re = Uy) (70) 
where 

U, is the free flow velocity beyond the influence of boundaries 

L is the length scale associated with the flow 


v is the kinematic viscosity of the fluid = 1.79 x 10-6 m-’s for water at 1 atmosphere 
and 0°C, 
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Whenever the Reynolds number exceeds a limiting value the flow becomes 
unstable and turbulence results. The instability is caused by the dominance of kinetic 
forces in the flow over viscous forces. For flow under ice, where relative current veloc 
tres can be of the order 0.1 to 0.2 m.s and length scales measured in thousands of metres, 
Reynolds numbers can easily exceed 2000, a typical limit for the onset of turbulent flow. 

In the presence of a rigid boundary, which in the first approximation can 
be considered flat and smooth, a turbulent boundary laver is formed. The thickness 6 
of this boundary laver is a function of the Reynolds number for the flow (White ,1986) 
such that 


é 0.16 


ee 7 
Ree 





(71) 


where 
x is the downstream distance from the begining of the boundary which is analagous 
to the lengtly Seale *eior ice sheetico.er 
If the boundary layer under the ice cover is fully turbulent it can be cone 
sidered identical to the mixed laver in terms of constant temperature distribution. Thus 
a tvpical value of 30 m, representative of arctic mixed layer depths, can be takenmi@nma 
the boundarv laver depth. It is accepted that processes other than turbulence in the 
horizontal flow contribute to mixing in the boundary layer under ice but the analogy 1s 


presented as a first approximation. From (70) and (71) 


Lo \tj6f 6 _ \76 + 
pose (£)"(a) : 


~ 3000m for Uy = 0.2ms7' and 6 = 30m 


This confirms that a rough length scale of the order of thousands of metres with a 
boundary layer depth of about 30 m is reasonable for large smooth ice sheets in wind 


driven motion. 
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b. Notse Mechanisms tn the Turbulent Boundary Layer. 

Vecchio & Wiley (1973) identified two mechanisms by which noise could 
be radiated from a turbulent boundary laver. The first is called Flow-Induced noise and 
involves acoustic dipoles induced in the boundary by turbulent pressure fluctuations. 
The second 1s called Flow-Excited noise and depends on the elastic behaviour of the 
boundary which, when excited by turbulent pressure fluctuations, radiates energy back 
into the fluid. 

The direct effect of the pressure fluctuations themselves, set up by turbulent 
motion in the boundary layer and often referred to as Psuedo-Sound or Quasi-Sound. 
could also be considered as a noise producing mechanism in its own right. Pseudo-sound 
1S a non-radiating phenomena. it’s direct effect only being felt when pressure fluctuations 
impinge directly on the active face of a hydrophone. It 1s closely related to self-noise but 
1s considered here to be purely a product of the turbulent flow and not dependent on 
hydrophone characteristics. 

Since the first two mechanisms mentioned above depend on the presence 
of pseudo-sound this phenomena will be examined first. both as a driving force for the 
other mechanisins and as a source of noise in itself. 

2. Pseudo-Sound. 
a. The Snooth Boundary Approximation. 

Skudrzvk and Haddle (1960) describe the length scales of turbulent motion 
in the boundary laver as varving between the thickness of the boundarv laver and the 
Kolmogorov microscale some 40 times smaller. Patches of turbulence are envisaged 
passing over a hydrophone or presumably any given point on the boundary at a rate 
eee ils 15 the convection velocity which is proportional to U, the free flow velocity such 
tat O,~0,SU,. The characteristic frequency /, of the resulting pseudo-sound pressure 


fluctuations is therefore given by 





U 
ines 0S —— a) 


where 


d is the separation between the patches of turbulence being convected past a given 
point. 


oT 


If d is approximately equal to the boundarv layer depth 6 
Up 
vip a (74) 


Thus, f,= 00007 Hz for U,=0.2 nis and 0 = 30mm 

Skudrzvk and Haddle (1960) found that the power spectrum of pseudo- 
sound pressure per unit frequency is virtually constant for frequencies up to f, and then 
decreases inversely proportional to a power m, at frequencies greater than f. Atf=f, 


the pressure due to pseudo-sound 1s given by 
Py = ().75 x 107°o.p2.U; 5* = (mn —1/m) Pa’ (7) 


where 
w% is the Kraichnan constant 
Gy 1s the free stream. clocity; 
o* 1s the displacement thickness of the boundary laver = 0/5 


p, 18 the water densitv 


So the power spectrum for pseudo-sound ts given by 


P(f) = 1.5 x 1077a" pid = (m — Ifn) Pa Hz fsh Ge 
and 
PU) = 1.5 x 10°82? ple 55 > (on — 1m) LN” PaHz7 f>fae 


o is related to the Revnolds number for the flow and hes between about 0.7 and 6, also 
Skudrzyk & Haddle (1960) give m= 3. Therefore, for water with p, = 1026 kg m° 


Pi) = 1.6U30°5 Pa*Hz7' sh GS 
and 
P(f) = 1.6U6075 ify Pa? Hz7! f>f aa 


ty 
Co 


Lotsch (1971) gives an alternate formula for the pseudo-sound power spec- 


trum in terms of the Strouhal number S 





P(f) = 3.16 x 10° p.6* UG + (2S)? Pa? Hz ($0) 
mere 

_ fe 

= 2 


n|o 


Ortakine 0* = andeas — 1026 ko mi 


fro ) —3/2 


Ss 





P(f) = 6.76U,( 1+ ( Pa°Hz ($1) 


Pressure fluctuations or pseudo-sound could affect a hydrophone located in the turbu- 
lent boundary laver and be recorded as noise. The estimated spectra. as given by (78) 
and (79) or (81), when compared to the measured spectrum from Makris and Dver 
(1986) show that pseudo-sound 1s probably insignificant as a noise source if the under 
sme O! thie ice 1s smooth (see Figure 33). Although the slope of the roll-off with in- 
creased frequency fits well with the measured data, the overall intensity is too small 


compared with the measured spectrum below | Hz. 
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Figure 33. Comparison of the calculated pseudo-sound specta (assuming no 





roughness) with a measured Arctic noise spectrum. 
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b. The Effect of Surface Roughness. 

The estimates above have assumed the ice sheet to be a smooth flat plate. 
This is a poor assumption as considerable underside roughness is known to exist 
(Wadhams. 1988). A commonly assumed value for the roughness scale under ice 1s 1 
m although higher values could be envisaged in some circumstances. 

The effect of surface roughness in boundary layer flow is to pierce the 
laminar sub-layer close to the surface of the boundary. Vortices are shed which pass into 
the turbulent region, increasing the turbulent motion there and generally shifting noise 
levels to higher frequencies. These effects begin to be important when the Reynolds 
number for the flow at the roughness scale exceeds a value of about 5 (Skudrzvk and 
Haddle, 1960), 1.e., when 
gas 


y 





ee fo.) 


where 


fees tie shear velocity, roughly representative of the velocity in the immediate vicin- 
itv of the Jaminar sub-laver, U* ~ 0.04U, (Skudrzvk and Haddle, 1960). 


and 


/iis the length scale for the roughness. 


mone — | in and U,=0.2 ms"! 


ao. -.- . 
= = 5000 (83) 
Therefore the rough underside of the ice can be expected to play a major role in the 
creation of turbulence. The characteristic frequency for pseudo-sound created by 
roughness 1s given by 


SL, 
Hatrouch) = h ($4) 


ad 


0. | ae ; 
Rather than the a given by (74) for the flat plate approximation. Thus taking. 
U,=0.2ms? and h=1 m, f,= 0.2 Hz is probably a more reasonable estimate for 


pseudo-sound under rough ice. 
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Equations (78), (79) and ($1) can now be rewritten for rough boundary 
conditions. The Skudrzvk and Haddle equations (78) and (79) become 


Pf) = 1.6)U2%75 Pa?’Hz7! IS focoum ae 
and 

ae ee 3 = 
P(f) = 1.6; US 076 Woe? Ba ia f > Socousn 


The Lotsch equation ($1) becomes 


ca 
Pf) = srs +( sah ae Pa*Hz (87) 


Step 





where 


y 1s a factor to take into account the increase in the intensity of the turbulent pressure 
fluctuations due to boundary roughness. 


The next problem is to estimate y. If the noise pressure 1s broadly propor- 
tional to the drag between the boundary and the fluid flowing over it, as suggested bv 
White (1979) y will be approximately equal to the ratio between the dimensionless drag 
coefficient in the rough condition and that in the smooth condition. White (1979) gives 
empirical equations for the dimensionless drag coefficient C, of a flat plate in both 
conditions and also presents the relationships graphically (reproduced as Figure 34). 


Chemoom 18 2 function of the Reynolds numberfor the flow:such that 


V7 


eC i 


0.03] 5 f Vol \aas7 
Crismooth) = = 0.031( 1 ) : (88) 


Cproughy 1S Independent of Reynolds number in the fully rough condition, depending in- 


stead on the parameter L,h (see Figure 34). 
Bercy = (1-89 + 1.62 log Lik)” (S9) 


L is the horizontal length scale for flow over the plate 


h is the roughness scale 


An estimate of the increase in pseudo-sound amplitude for a rough surface will therefore 


evctyen by the ratio 


ee : 


(lL) = = (90) 
Crismooth) 0.031(1.89 + 1.62 log Lik)" 


Note that this is only valid for fully rough conditions, 1.e., the transition y — unity for 


fi — 01s not described. 
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Figure 34. Dimensionless drag coefficients for rough and smooth flat plates (from 
White, 1979). 


Figure 35 shows the calculated pseudo-sound spectra from equations ($5) 
(S86) and (87) for a roughness height h =1m, a boundarv layer depth 6 = 30 m and a 
free flow velocitv U, = 0.2 ms-!. As with the smooth plate approximation, the general 
shape of the curves fit well with the measured spectrum below | Hz. The characteristic 
frequency f, , however, is a little too high suggesting that the effective length scale may 
In reality lie somewhere between 6 and h. Despite this the fit looks pronusing given the 
approximate nature of the theory used, and the fact that the input pararameters are 
crude estimations. This indicates that this mechanism is probably responsible for a large 
proportion of the noise at these low frequencies. It is worth noting at this point that 
Makris and Dyer (1986), in a note accompanying thelr composite ambient noise spec- 
trum, hypothesise that the noise below | Hz may be due to non-linear surface wave noise 
or psecudo-sound in the oceanic boundary laver. No attempt is made to fit a curve ex- 
peeeevo the measured data, as none of the free variables, U, . 2 or A are known for the 
data set used to produce the composite spectrum. 

In conclusion, the indications are that pscudo-sound may well be an 1m- 
portant source of noise for hvdrophones within the turbulent boundary laver under 
how@eueice s1.c.. positioned io more than a few tens of metres below the ice. The fre- 
M@eeneies involved are. however. generally below | I1z, a region which 1s not, at present, 


of general importance. 
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Figure 35. Comparison of the calculated pseudo-sound specta (roughness effects in- 
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cluded) with a measured Arctic noise spectrum. 
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3. Flow-Induced Noise. 

This is a mechanism by which the effects of turbulent pressure fluctuations are 
radiated bevond the region of turbulence to the far field. The effect 1s due to two phe- 
nomena. The first occurs within the volume of the turbulent fluid and was described by 
Lighthill (1954) as being due to quadrapole sources. The second occurs at the interface 
with a sohd boundary and is described as being dipole in nature (Lotsch, 1971). Acoustic 
power due to the first of these effects 1s proportional to the eighth power of the Mach 
number while the power due to the second is proportional to the Mach number squared 
(Vecchio and Wiley, 1973). In the under-ice environment Mach numbers are very small, 
typically about 10-*. The dipole sources therefore dominate but the total radiated power 
is considerably smalier than that due to the turbulent pressure fluctuations themselves. 
Spherical spreading further reduces sound intensities due to this radiative mechanism. 
Flow-induced noise is not, therefore, considered to be an important contributor to 
under-ice noise. 

4. Flow-Excited noise. 

Vecchio and Wilev (1973) identified this mechanism as being another method 
by which noise can be radiated from a turbulent boundary. The mechanism depends on 
the presence of turbulent pressure fluctuations along a flexible, elastic boundary. 
Vibrational modes are excited in the boundary which then radiates sound back into the 
body of the fluid. Floating ice can certainly be considered as a flexible, elastic boundary; 
also a preceding section has shown that low frequency pressure fluctuations are created 
i 110w Over its roush surface. [here exists a possibility. therefore, that this effect might 
be responsible for noise beyond the turbulent boundary laver nmmediately under the ice. 
The physics involved are, however, extremely complex. For this reason no further in- 
vestigation has been undertaken. Detailed analvsis of this process is Warranted and mav 


lead to the conclusion that this is a mechanism of some significance. 
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5. Resonant Cavities. 

Water flow past cavities in the under side of an ice sheet is a potential source 
of noise which must be exanuned in order to establish its possible inportance in the 
overall creation of ambient noise. Open cavities or orifices in solid boundartes are known 
to resonate under certain conditions. The onset of resonance depends on the fluid flow 
velocity past the cavity, the fluid characteristics and the shape and size of the cavity it- 
self. There is water flow past cavities in the under-surface of ice sheets or the submerged 
surfaces of pressure ridge keels whenever a relative current exists. The object of this 1n- 
vestigation is to determine whether this flow is likely to be sufficient to cause noise from 
such cavities. 

a. The Helinholtz Resonator Model. 

The Helmholtz resonator is a rigid walled cavity whose resonant properties 
are known. Theory pertaining particularly to the frequency of resonance as a function 
of cavity dimensions is used here to model under-ice cavities. Figure 36 shows cross- 
sections of the classical Helmholtz resonator and an idealized simple cylindrical cavity 


such as nught exist in the under surface of an ice sheet or in a pressure ridge keel. 


Classical Helmholtz 
Resonator 


de 





Figure 36. The classical Helmholtz resonantor and a simple cavity. 
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Vis the volume of the inner cavity which provides the mass element of the 
resonator. S: is the area of the circular opening which acts as a simple acoustical source 
and the fluid in the neck, length L. provides the sziffvess or spring component. From 


insier et al. (1952) the frequency of resonance 1s given bv 


a ~ C( 757)" (91) 


C is the sound speed of the fluid 


L’ 1s the effective length of the neck taking radiation -mass loading into account. For 
practical purposes L’> L+0.8A 


A is the diameter of the opening 


It is assumed that the acoustic wavelength 4 of the sound resulting from resonance 1s 


Ieee Compared with the length scales L, S'? and #3. For the idealized ice cavity 


L=B-+ (92) 

oa <( = ) (93) 

and 

ceed x( oi (94) 
hak 

Therefore 


= C 1S i 
tein = of 


fo taking C = 1500 ms~ for sea water 
Pee 5810 (104B — 34°)! 
and 


fo =~ 925 (10AB - 34°)? Hz (95) 
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The constraint on / relative to the cavity dimensions dictates that A ~ B for (95) to be 


valid. Figure 37 shows resonant frequencies for such simple cavities. 


0.00 0.5 
Dimension B (metres) 


Figure 37. Resonant frequencies (Hz) for simple cavities. 





b. Excitation of Resonance. 
The next concern ts whether the modest flow velocities present under ice are 


suflicient to excite resonance in the cavities there. From Blake (1984) the relationship 
between the free-flow velocity across the entrance to a cavity and ff the frequency of the 


pressure fluctuations caused by the flow is given by 


ee 
a fe 1/4) n=1,2,3,.... (96) 


U, is the free flow fluid velocity 
C.is the hydrodynanuc phase velocity across the opening 


A is the diameter of the opening 


ee CG 
For turbulent boundary layer excitations 7 = Ono SO 
0 


0.33U 
f= eras (a— 1/4) n= 1,2,3,... (97) 
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Figure 38 shows excitation frequencies /; for cavities up to one metre in diameter. A 
nonunal value of U, = 0.2 ms-' has been taken to represent a typical maximum relative 


current velocity under ice. 
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Figure 38. Excitation frequencies for siniple cavities (Hz) in a flow velocity of 0.2 





m/s. 


By comparing Figures 37 and 38 it can be seen that the frequencies at which 
under-ice cavities could be excited by a typical flow velocity are far lower than their 
resonant frequencies. 

Thus, it 1s concluded that, within the limits of the approximations inade, 
resonant cavities are unlikely to be important in the generation of ambient noise under 


ice. 
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V. SUMMARY AND CONCLUSIONS 


The importance of being able to accurately predict ambient noise levels in Arctic 
waters was discussed in Chapter I. In order to attain that abilitv research must be 
undertaken into every aspect of the Arctic ambient noise problem. The reviews of 
Chapter I] and the statements made by Pritehard (1988) illustrate the great number of 
different subsidiary problems which must be understood and overcome before a reason- 
able solution can be reached. The large number of possible mechanisms involved in noise 
generation and their extreme complexity make this a verv daunting task. In comparison 
to the problem of modelling open ocean noise, the Arctic problem 1s complicated by the 
different interactive processes between air, 1ce and sea. In addition, the inaccessabilitv 
of the Polar regrons leads to a dearth of meaningful data. Progress 1s, however, being 
made. Wenz (1962) commented on the inadequacy of data from open ocean sites at that 
time, but was able to describe and catalogue the nolse-generating processes to good ef- 
fect. A stmular situation exists now with regard to noise produced in 1ce-covered waters. 

Pritchard (1988) categorized Arctic noise mechanisms into three types bY foram 
function (as described in Chapter I). This work addresses the mechanisms constituting 
the third category, 1.e., those mechanisms dependent on the shearing of the boundary 
laver of the ocean under the ice. A scheme is proposed in Chapter II! by which this 
category could be further sub-divided into Flow and Flow’ Mechanical mechanisms. The 
following is a svnopsis of the results of investigations into this sub-set of potential notse 


mechanisms. 


A. FLOW/MECHANICAL MECHANISMS. 
I. Cracking Noise Due to Wind/Current-Induced Moments. 

From statistical work done with observed data, Makris and Dver (1986) found 
that ambient noise correlated well with the stress moment acting about an ice sheet's 
central horizontal plane. In this study, bending moments induced tn the ice sheet were 
investigated by considering the effects of wind and current on an idealized pressure ridge 
feature. Ridges were envisaged as being in one of two conditions, the young unconsol- 
idated ridge and the more mature consolidated feature. In the case of the first it was 
found that the moments caused by wind and current are insignificant when compared 
to the bending moments caused by isostatic loading. The relationship between the 


bending moments present due to isostatic loading and the critical bending moment re- 


quired for cracking in the ice sheet was identified and unstable combinations of ridge 
height and ice thickness were found. There appeared to be a remote possibility that the 
wind and current-induced bending moments might act as a precipitative mechanism for 
cracking actually resulting from isostatic loading. It is felt, however, that this 1s not a 
likely explaination for the correlation observed by Makris and Dyer (1986). The study 
of the consolidated ridge revealed even less likelihood of wind’current-induced cracking 
through induced bending moments. The assumed rigidity of the system resulted in 
bending moments that were far below the critical value required for initial cracking. 
2. Ice Fragment Bumping 

Beyond some observations that ice debris does seem to move about in current- 
driven motion causing audible sounds, very little appears to have been written on this 
subject. The concept is attractive in that similar ice-ice interactions are known to cause 
a great deal of noise at the ice edge. Given sufficient amounts of freely moving fragments 
and a moderate relative current, a verv noisy result could be envisioned. Although the 
detailed physics of this processes were not investigated, there are indications that the 
resulting sound would be spread over a wide frequency range, possiblv varving from tens 
of Hz to the low kHz. This is a pronusing mechanism warranting further investigation. 
The problem of modelling this effect in practical terms, however, does present many 
difficulties. The most accute would probably be the determination of the density of free 


ice fragments in any given region. 


B. FLOW MECHANISMS 
I. Noise From Turbulence 

This is the process by which noise is produced as ice moves in Water or when 
water flows under static ice. It 1s probably the process envisioned by Lewis and Denner 
(1988) when thev wrote of noise caused by ice “rushing through the water.” Investi- 
gations of the turbulent pressure fluctuations which could be expected for typical ice 
roughnesses and flow regimes revealed a predominance of very low frequency activity. 
When compared to a measured Arctic spectrum, the predicted turbulent pressure or 
pseudo-sound spectra have the correct characteristics to explain low frequency noise. 
Thus it 1s suspected that pseudo-sound in the boundaryv layer is significant at frequencies 
below | Hz but only for receivers actually within the region of turbulence. 

Two mechanisms, Flow-Induced and Flow-Excited noise, by which the effects 
of turbulent pressure fluctuations might be radiated were mentioned briefly. The first 


appears to be too inefficient to radiate significant sound energy. The second may, how- 


? 


ie 


ever, be of some importance. Both phenomena need further, more detailed, consider- 
ation. 
2. Resonant Cavities 

The idea that cavities in the under surface of the ice night resonate due to water 
flow comes from a consideration of simular features in air. A comparison with the clas- 
sical Helmholtz resonator concept revealed the likely frequencies of resonance for vari- 
ous cavity dimensions. By establishing the frequencies that tvpical flow velocities might 
excite, however, it was found that such flows are almost certainly too slow to cause 
resonance. Thus this mechanism is considered to be an unlikely candidate for noise 


production. 
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APPENDIX BEAMS ON ELASTIC FOUNDATIONS. 


This appendix describes the theoretical behaviour of infinite and semi-infinite beams 


On elastic foundations as given bv Hetenvi (1946). This theory is used to model floating 


ice because the ice 1s flexible to a degree and since buoyancy and gravity make the sea 


act like an elastic foundation. 


A. INFINITE BEAMIS. 


The spatial distribution of displacement, v(x), slope @(x), bending moment J/(x) and 


shear force O(x) for an infinite beam, on an elastic foundation, subjected to a point 


foading P and moment Jf (see Figure 39) are given by:- 











P} Ngee: 
ie) = 7 Bye Se F Bix 
Pji* We 
Ox) =--*- B+ C,, 
jE M/ 
Be —- Ct Oae 





where 


Ais the foundation modulus =(p,.— p)g for floating ice 


a; Ne oa oe 
iis the Characteristic length of the system = Ge m| 


Eis the Youngs modulus for the material of the beam per unit width 


3 
l 
I 1s the moment of inertia for the beam =-—— 


J 


ris the thickness of the bean. 


~~] 
cn 


(47) 


A), = e~*( cos Ax + sin Ax) 


—J/xX . 
By.=e ~ sindx 
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Ch, =e cos Ax — sin /x) 


Ax 


D,, =e ~ cos Ax 


(48) 
(A9) 


(A 10) 


(A11) 





Figure 39. The infinite beam with a point load and bending moment. 


B. INFINITE BEAMS WITH TRIANGULARLY DISTRIBUTED LOADINGS. 


For a triangularly distributed loading (see Figure 40) 


q 
u(x) = ae (Cy) px) — Cox — 24D), + 4H(E — x)) 


Jo 
WO) = Gage (Crienet = Cox = 22€D,,) 


0 
OE ay eas Ce ee ee 
O(x) = -— (By ,_y + By — 2LCy,) 
ape? 
where 


gis the load at x =0 


f is the base length of triangular loading 
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Figure 40. Infinite beam with triangularly distributed loading. 


C. SEMI-INFINITE BEAMS WITH TRIANGULARLY DISTRIBUTED 
LOADINGS. 

If the infinite beam equations, (A1) to (A4), are to be used to describe a semi-infinite 
beam with a triangular loading, end conditioning forces must be applied in order to 
create the effect of a free end. Thus, if a free end is to be created at x = 0, a force P’ and 
moment J/’ must be applied to exactly cancel Q, and Af, The result will then be 
‘1(0) = 0 and QO(0) = 0, the conditions for a free end. From (A14) and (AI5) at x =0 


do 
a ee (A16) 
0 giap wnat ) 
Qy = ——2~ (By, — 20) (417) 
aize 
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Substituting (A16) and (A17) into (A3) and (A4), with x =0 





xO) = ay re ; 2p ( wali=9 
— P’ AM’? q 
90) ==--=+- o> (Bre — 44) =0 


Solving (A18) and (A19) for P’ and Jf’ 








qo : 
P= (al os eae a 
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(419) 


(420) 


The distribution of displacement on the semi-infinite beam with a triangular loading 


can now be found by superimposing the effect of the load from (A12). together with the 


effects of Panduit. = 01s described a acum): 


ae 





ene a 

(x) = Mia rare Fan Cr = 2A¢D), a8 a" es) ne Aj + 
Io WDD M' 

WO) =e (Gal — Ox Oe 


Substituting for P’ and J/’ and simplifving 


do 
i. dkAo (Capex) FAA —x) + By Ane — Dj (Age + 2832) 
ACS) me = 5 Creel Oe ee) 


ay 
And since J/(x) = —- Ela (from Fleteny1) 


q 
M(x) =- ur (Aji eey| ~ Die Ee 


— By, OS x<¢(A22) 
x > ¢(A23) 

O<x< (A224) 
x>C(ee 

x > 0(A26) 


D. SEMI-INFINITE BEAMS WITH POINT LOADINGS AND BENDING 
NIOMENTS. 

Developments similar to those describing the behaviour of a semi-infinite beam with 
a triangularly distributed loading are used in the case of a point loading and bending 


moment at the free end (see Figure 41). 





Figure 41. Semmi-infinite beam with a point loading and bending moment. 


Equations for the displacement, slope, bending moment and shear force distibutions 


are given below. 





2P) SVE 
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